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Abstract 

In the present work we study an anisotropic layered superconducting film 
of finite thickness. The film surfaces are considered parallel to the be face of 
the crystal. The vortex lines are oriented perpendicular to the film surfaces 
and parallel to the superconducting planes. We calculate the local field and 
the London free energy for this geometry. Our calculation is a generalization 
of previous works where the sample is taken as a semi-infinite superconduc- 
tor. As an application of this theory we investigate the flux spreading at the 

superconducting surface. 

PACS numbers: 74.60.Ec, 74.60.Ge 
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Scanning superconducting quantum interference device (SQUID) microscope has been 
used to image interlayer Josephson vortices trapped between the planes of layered super- 
conductors. This technique has been used to measure the out-plane London penetration 
depth that gives the distance over which the interlayer current j c changes as a function of 
in-plane coordinates.0H These measurements have been important to test the interlayer tun- 
neling model as a candidate to explain the mechanism of superconductivity for the high-T c 
superconductors.ifl 

Recently, Kirtley, Kogan, Clem, and Moleri have found expressions for the local magnetic 
field emerging from a superconductor with the vortex lines parallel to the planes, and normal 
to a crystal face. Their geometry consists of a semi-infinite anisotropic superconductor. 
Furthermore, they have used these expressions to fit the experimental data at the surface in 
order to obtain an estimate of the value of the out-plane penetration depth A c . They have 
shown that, neglecting the vortex spreading at the surface may overestimate A c as much as 
30%. 

In the present paper we extend the work of Ref. [5] to an anisotropic layered supercon- 
ducting film of finite thickness and of infinite extent in the be face of the crystal. We will 
show that, if the thickness of the film is of order or smaller than A c , the magnetic field 
distribution is even more affected by flux spreading. 

Let us first formulate the problem to be solved. The geometry we consider is illustrated 
in Fig. |l} We suppose that the vortex line is perpendicular to the film. We will calculate 
the local field inside the film using the London equation. For this geometry this equation is 
given by 

V x [V Xh] +h = z$ 5(r) , (1) 

where V is the London (tensor) penetration depth. This tensor is diagonal and its compo- 
nents are given by A xx = A yy = A^, A zz = A^ b ; here X a b and A c are the in- and out-plane 
penetration depth respectively; $o is the quantum flux. The film is anisotropic along the c 
direction. 
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Outside the sample, the local field satisfies the equation 



V 2 h = . 



(2) 



Although we will consider the case of a single vortex, the generalization to the case of N 
vortices is straightforward. To proceed is more convenient to Fourier transform Eqs. ([!]) and 
(@). For \z\ < d/2, using the Maxwell equation V • h = 0, we obtain a set of three coupled 
differential equations for the two dimensional Fourier transform of the local magnetic field 
h(k,z) =/ c/ 2 re- lk - r h(r,2) , 

d 2 ~ 



1 -t- A ah K A ab 



h T = 



(3) 



1 + \ ah k - A 



c dz 2 



h y + (A 2 - \ 2 ab )ik 



' dz 



0, 



(4) 



1 + A ahkl + A 2 A; 2 , — A 



a 2 

'dz 2 



K + (A 2 - X 2 Jik y ^ 



For | z | > d/2 one has 



dz 2 



k 2 )h = 



(5) 



(6) 



At the vacuum-superconductor interfaces z = ±d/2 the field components are continuous 
and the component of the current perpendicular to both film surfaces vanishes. One has 



h<(k, -d/2) = h m (k, -d/2) 



(7) 



h m (k,rf/2) = h>(k,ci/2) 



z ■ [D z (k) x h m ] z=±d/2 = , 



D 2 (k)-h = 



(8) 
(9) 
(10) 
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where the operator D z (k) = ik + zjp The subscripts (<, >) stand for below the surface 
z = —d/2 and above the surface z = d/2, respectively, whereas the subscript m is meant for 
the field inside the sample. 

We start by solving first Eq. @. The solution which satisfies the boundary condition of 
Eq. ([10]) takes the form 



h>(k, z) = (-ik + zk)<p(k)e- k{z - d/2) 
h<(k,z) = (ik + zk)<p(k)e k{z+d/2) , 



(11) 
(12) 



where <f{k) is a scalar function which will be determined by using the boundary condition 
either of Eq. (g) or <g). 

Eq. (g) can also be easily solved. We have 



h m>x (k,z) = W 1 e az + W 2 e 



(13) 



where 



a 



1 + Aj t P 



2 

ab 



(14) 



and the W's are two constants to be determined by using the boundary conditions. 

The other two components of the local field can be determined by decoupling Eqs. (f§) 
and (|j). This can be done by calculating the determinant of the matrix formed by the 
coefficients of Eqs. (f|) and Eq. (^). This yields the following equation for h m ^ 



,1x2.2 x2 9 2 
1 -+- A ab K A ab 



1,\2 1.2 , \2 h 2 x2 ® 2 
1 "r A ab^x "r A c^y c Q z 2 



(15) 



The solution for this equation is given by 



h y>m {k, z) = W 3 e az + W 4 e~ az + W^ z + W 6 e 



(16) 



where the W's are constants to be determined by using the boundary conditions and 
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\ 



1 + A ab^x + A c^y 

X 2 . 



(17) 



The solution for h m>z can be found by inserting Eq. (|T6| ) back into Eq. @ or (|jJ).One has 

(18) 



A c 7 iK y 7 



The determination of the constants Wi is very cumbersome and we omit it here. We just 
present the main steps of the complete solution. First of all, we use the Maxwell equation 
(flOf). This allows us to write W 3 and W4 in terms of W\ and W 2 - Secondly we use the 
boundary condition of Eq. (||) in both faces of the film. This leads us to the solution of W5 
and Wq in terms of W\ and W%. Then, we are left only with three constants to determine, 
namely, W\, W 2 , and ip. Thirdly, we use the continuity of the local field at the film surfaces 
[either Eq. (0) or @; both of them yields the same solution to these constants]. One obtains, 



Wo = ik n 



W 3 

w 5 

where 



2 sinh (f^) 
2 sinh (^f) ' 



1 + A ab^x 
1 + A ab^x 



2 ■ 



k x (l + \ 2 ab kl) sinhfrd) 



k x (l + \ 2 ab k 2 x ) smh(jd) 
$0 



W\ sinh 
W\ sinh 



(7 + «) 2 



(7 - «) 



d 



+ W 2 sinh 
+ W 2 sinh 



(7 - a) - 



(a + 7) 



d 



A c 2 7 2 



A(k), 



(19) 
(20) 

(21) 
(22) 
(23) 
(24) 
(25) 



A(k) 



k + 



A^«coth(f) + f coth( 



1 + Kbkl 



(26) 



Finally, upon substituting Eqs. (p~9| - |25D into Eqs. ([13|), 
magnetic field inside the film 



^m,a:(k, Z) 



-ik x if(k) 



sinh(az) 
sinh (ef) 



and (1181), we find for the local 



(27) 
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^m,2/(k, z) — iky 



vp(k) 



1 + ^ab^x 



, 9>9 sinh(az) sinh(72:) 
XLk 2 



ab x 



+ 



sinh (f) sinh 



Ah 2 



A 2 fr 2 



,,,9 cosh(a^) A;, 2 cosh(7z) 



+ 



(28) 
(29) 



sinh(f) 7 sinh (f) 
We would like to point out that these results could not be obtained from those of Ref. || 
without solving the problem. In fact, the solution of the London equation for a supercon- 
ducting film is different and more difficult than for a semi-infinite superconductor. 

Let us turn our discussion to the calculation of the London free energy. The energy of 
the vortex system is given by F = Fy + Fg, where Fy is the field energy in the vacuum and 
F$ is the energy inside the superconductor. One has 

1 r d 2 k { r° . .. .. ... r- d / 2 



Fy 



57T 



57T 



(2tt) 2 

d 2 k 



d/2 
d/2 

-d/2 



dz |/i>(k, z)\ 2 
dz {|/i m (k,z)| 



r- a/2 

/ dz\h < (k, z) 

J —oo 



+ [D,(k) x h OT (k,^)] ■ V • [D 2 (-k) x h m (-k,z)] 



(30) 



(31) 



By substituting the appropriate expressions of the local magnetic field inside Eqs. (|30] 
and (|3l|) , after a length algebra, we obtain 



$5 



d 2 k 



d + 2 



A(k) 



Ah 



(32) 



8tt J (2tt) 2 A 2 7 2 

The free energy can be generalized to an ensemble of N interacting vortex lines upon 
multiplying the integrand of Eq. (|32|) by ^(k)) 2 where the structure factor is given by 



S(k) = £ e 



tk-R, 



(33) 



Here R« is the position of the i-vortex line. Note that this extended result should be 
valid for an ensemble of distorted vortices, that is, the positions of the vortices do not 
necessarily correspond to the equilibrium configuration. The first term inside Eq. ( |32| ) 
represents the interaction energy of the vortex lines as if the surfaces were absent. The 
second term represents the surface energy associated to the magnetic energy of the stray 
field at the superconductor- vacuum interface. Notice that for k small (large r), 7 2 ~ 1/A 2 , 
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and A(k) ~ 1/k. Thus, the surface energy goes as <£>Q/87r 2 r. Consequently, the interaction 
on the surface depends neither on the film thickness nor on the anisotropy. This is the 
Pearl result for vortices emerging from a semi-infinite isotropic supercondutor.il Another 
interesting particular case of Eq. (^) is the limit of a very thin film d — > 0, and k small. In 
this limit, from Eq. ( PBj ) it is straightforward to show that A(k) = 1/ (k + 2(\ 2 ab k 2 .+ \ 2 k 2 )/d). 
Therefore, from Eq. ( |32"|) we obtain 

f dk 2 2nd 
J (27p fcA- 1 + (k 2 + Tk 2 ) ' ( ' 

where Eq = ($o/4vrA a fc) 2 , A = 2X 2 b /d, and T = \ 2 /X 2 ab is the anisotropy parameter. This is 
precisely the energy of a single vortex in very thin film first obtained by Pearl.Bi 

Now we will turn our attention to the streamlines of the integrated field over x. The 
distribution of magnetic field emerging on the surface can be probed with a SQUID pickup 
loop. If the SQUID probe is oriented in the xy plane, the total magnetic flux will be nearly 
equal to the pickup loop size timesi 

/oo /*oo A lr 

h z (x, y,z)dx= — * h z (0, k y , z)e ik yy , (35) 
-oo J —oo Z7T 

whereas, if the SQUID probe is oriented along the xz plane, the total magnetic flux is 
measured through the pickup loop size times 

poo poo /jU 

H y (y, z) = I h y (x, y,z)dx= I h y (0, k y , z)^ . (36) 



oo J — oo 



In order to compare our results with the results of Ref. ||, we will replace the vacuum- 
superconductor surfaces at z = and z = —d. This can be done through the translation 
z -> z + d/2. From Eqs. (0), (0), and (H^H) we obtain, 



h> (k, z) = (-ik + zk)<p{k)e- kz , (37) 
h< (k, z) = (ik + zk)<p(k)e k(z+d) , (38) 



. ik ^ I X 2 fc 2 



sinh 



a 



sinh 



<^(k) 



A 2 7 2 1 + XLkl 



X 2 ab k 2 x a- 



sinh (?) 
cosh 



+ 
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a 



(*+f) 



sinh (?) 



sinh (?) 



+ 



A; 2 COsh 

2/ 
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M) 



7 



sinh (?) 



(40) 
(41) 



The substitution of the appropriate expressions into Eqs. ( |35| ) and (|36| ) yields for the z 
component of the 7i field, 

$ n r 00 , cos (?/' sinh m) e ~ z ' sinhu 



7rA c JO 



cosh u + sinh u coth (^- cosh 



it 



(42) 



$ (TV 



TrX r 



x 



2 6 



tanhw 



cos(?/sinh w) 


cosh 


l( Z ' + A) 


coshw 


1 


cosh m + sinh u coth (^- cosh w 


sinh (^- coshw 


J 



00 , cos(VsinhM)e (2 ' +d/Ac)sinh?t 
riu K - '- 

cosh u + sinh u coth f cosh u 



7rA c jo 

where = y/X c and z' = z/X c . 

The y component takes the form 

$n Z" 00 , sin (2/ sinh w) e" z ' smh 



(43) 
(44) 



ni{ y ,z) 



_jo_ 

7l\ c Jo 
$ 



cosh u + sinh w coth cosh « 



7rA c jo 
$ 



du 



°° , sin(!/sinhw)e (z ' +d/Ac)sinhu 
dtt- 



sin(?/ sinhw) 


sinh 




+ 2^-) coshw 


cosh u + sinh m coth ( cosh w 


sinh 


(_£ coshM ) 



(45) 

(46) 
(47) 



n\ c Jo cosh u + sinh u coth (2^- cosh 
Note that in the limit of d — > 00, our results are exactly the same as those of Ref. |5|. 
The results for the 7i field presented above should be useful to interpret the experimental 
data obtained by using scanning SQUID micorscopy. Unfortunately, the experiments have 
been performed in samples of large thickness.0! This renders the test of the theory impracti- 
cable. In fact, vortices have been magnetically imaged in films, but for a different geometry, 
that is, the superconducing planes are taken parallel to the surfaces of the film and the 
vortex lines are considered perpendicular to the film surfaces.^ In this case, we can extract 
the in-plane penetration depth A a t, rather than A c , from the fitting of the experimental data. 
So, we will restrict our analysis only to the theoretical expressions. 
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Fig. |2] shows the streamlines of the T~C(y, z) field for a single interlayer vortex centered 
at x — 0, y = 0. The streamlines were generated as sketched in Ref. ||. We used various 
values of the film thickness. Note that as the thickness of the film grows, the flux spreading 
is important only near the surface, whereas deep inside the thinner film the streamlines are 
still very distorted, except those close to the center of the vortex. 

To see how important the flux spreading inside a superconducting film is, we calculated 
numerically ir\ c TC z (y, z)/Qq as function of y/X c for three different values of d at z = 0. As 
can be seen from Fig. [|, the full width at half maximum of the flux contour is 1.87A C for 
the case d = 5A C , while it is 1.65A C for d = A c . Thus, if the flux spreading inside the film is 
not taken into account, the value of A c could be underestimated by 10%. This error grows 
as the film thickness decreases. 

Finally, we would like to point out that the present results agree with their isotropic 
counterpart. If we set \ a b = A c = A in Eq. (p5|) and fl26|) , we obtain the same result as in 
Ref. H Apparently, our results are different of those found in Ref. |10|, but they show very 
similar streamlines. 

In summary, we have calculated the field distribution of a single vortex inside and outside 
a layered superconducting film of arbitrary thickness. We also calculated the London free 
energy of an ensemble of vortices. From the expression for the energy one can recover the 
interaction potential between vortices for a very thin film0i and the vortices emerging from 
a semi-infinite superconductor .@i In addition, we have shown that flux spreading inside a 
superconducting film of order or smaller than A c affects substantially the full width at half 
maximum of the flux contour. 
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FIGURES 

FIG. 1. Geometry of the film used in this work. The vortex lines are oriented perpendicular to 
the be face of the crystal. 

FIG. 2. Streamline mapping of the integrated field H, for an anisotropic superconducting film. 
The spacing of the streamlines is proportional to (dHz/dy)^ 1 at z = — (d + 0.5A C ). (a) shows the 
streamlines for a film of thickness d = A c , (b) for d = 2.5A C , and (c) for d = 5A C . 

FIG. 3. The z component of the integrated field as a function of y at z = for three different 
values of the film thickness: (a) d = 5A C (dashed line); (b) d = 1.5A C (dot-dashed line); (c) d = A c 
(continuous line) 



11 




a, z 



d 



Fig. 1/Sardella 



12 
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